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ABSTRACT 

In this paper we present a detailed study of the peculiar early-type galaxy 
NGC1947. The main goal of this work is to constrain the dynamical status 
and the formation history of NGC1947 by comparing the observed properties 
with the predictions derived from different galaxy formation scenarios. To 
this aim, we derived the photometric and kinematical properties of NGC1947. 
Due to the presence of an extended dust-lane, which crosses the galaxy center 
along the photometric minor axis, we used near-infrared images (J and K 
bands) to derive an accurate analysis of the stellar light distribution. Optical 
images (in the V and R bands) are used to derive the color profiles and 
color maps to study the structure of the dust-lane. The observed kinematics 
confirm the presence of two components with decoupled angular momentum: 
gas and dust rotate along the minor axis, while the rotation velocities of the 
stars are observed along the major axis. The complex structure observed in 
NGC1947 support the hypothesis that some kind of interactions happened in 
the evolution of this object. We analyzed two alternatives: a merging process 
and an accretion event. We discussed how the observed properties strongly 
suggest that the decoupled ring of gas and dust have been accreted from 
outside. 

Key words: galaxies: individual - galaxies: formation - galaxies: evolution 
- galaxies: interactions. 
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1 INTRODUCTION 

One of the major open issues in modern cosmology is to understand how galaxies formed 
and evolved. It is likely that the formation of galaxies is dominated by two processes: the 
assembly of luminous and dark matter through accretion and mergers, and the conversion 
of baryonic and non-baryonic matter into stars. The Cold D ark Matter (CD M) scenario for 
galaxy formation is based on the hierarchical mass assembly (jCole et al.ll2000l ) which predicts 
that the observed galaxies and their dark halo were formed through a repeated merging 
process of smaller systems. This is the reason why the study of galaxy interactions and 
merging has received an ever increasing attention in recent years, both on the observational 
and theoretical sides. 

It has been shown that gravitational interactions and mergers affect the morphology and 



dynamics of galaxies from the Local Group to the high-redshift universe (e.g. 



2003 



Bundy et al. 



2004 



Lin et al. 



2004 



Kartaltepe et al. 



2007 



Conselice 



De Propris et al. 



2007 



et al 



Lotzeta 



20081 ). These observational results support the scenario in which the merging of two disk 



galaxies produces spheroidal merger remnants with physical characteristics, such as density 



profiles, m ean velocity dispersion and surface brightness, quite similar to t 



rose observed for 



E and SOs (iToomre fc Toomrdll972l ; iBarnes fc Hernquist 



1992 



Bekki 



1998a 



Bournaud et al 



20051 ). However, other types of interactions, such as smooth accretion and stripping seem 
to be equally important in the growth of galaxies and the relative weights of the various 
processes depends on the environ ment, which driye man y morphological features observed in 



galaxies, such as bars, s hells (see 



rotating rings/disks (see 



Galletta 



Malin fc Carter 



19801 ) and highly inclined or even counter- 



19961 for a review). Concerning the latter morphology, dif- 



ferent degrees of evolution of these peculiar rings/disks are seen: some cases the decoupled 



matter is compose c 



(llodice et al. 



2002 



by stars and gas, and regions of recent sta r formation are also observed 



Iodice et al. 



2002: 



a par ticularly smooth structure (llodice et al 



Reshetnikov et al. 



20051 ): other cases sh ow rings with 



and 



20041 ) , or formed only by gas ( 



20011). In general, t 



Galletta et al 



e.g. 



rese p eculiar galaxies have very large amount of gas (IRichter et al 



Bettoni et al 



1994 



Arnaboldi et al. 



1997) an d dust which lies in the same plane of t 
' II 1 I 



1997 



Bettoni et al.ll200ll : Ivan Driel et ai- 



r e decoupled ring / disk 



2002 



Cox et al 



20061 ). The 



decoupling of the angular momentum between the two components cannot be explained by 
the collapse of a single protogalactic cloud, thus a "second event" must have happened in 



* E-mail: spavone@na.infn.it (MS); iodice@na.astro.it (EI) 



Revisiting the formation history of the minor-axis dust lane galaxy NGC1947 3 

the formation history of these systems. Given all the above, such peculiar galaxies turn out 
to be the ideal laboratory where to derive detailed information about many of the physical 
processes controlling galaxy interactions and merging. 

Several theoretical works based on numerical simulations, have tested two plausible sce- 
narios for the origin of the decou pled matter: either gas accretion from infall processes, or 



merg i ng between gas-rich galaxies (IThakar fc Rvden 



1997 



Bckki 



1998b 



Bournaud fc Combes 



coupled material forms a ring ( 



Quinn 



200.4 



1991 



1996 



Thakar fc Rvde 



Bournaud et al 



„ 



1998 



Reshetnikov fc Sotnikove 



2005]). In both cases, the de- 



Hernquist fc Weill Il993l ) which after a while 



settles into one of the principa l planes of the gravitational potential associated with the 



host galaxy (IHeisler et al. 



1982 



Bertola et al 



1991 



Bournaud fc Combes 



2003|). Multiple 



accretion events may also account for more complicated objec ts, like NGC2685 where two 



orthogonal rings are obser ved around the main galaxy body (ISchinnerer fc Scoville 



Bournaud fc Combes 



2002 : 



20031 ). 



The unequal-mass merging of two disk galaxies is the other way to build a galaxy with an 
highly inclined ring/ disk: the morphology and kinematics of the merger remnants depends 
on the initial or bital parameters and the initial mass ratio of the two galaxies involved in 



the process (e.g. 



Bournaud fc Combes 



2003 



Bournaud et al 



20051 ). 



External gas can also be accreted from cosmi c web filaments w ith high inclined angular 



momentum relativ e to the main disk plane (IDave et al. 



2001 



Semelin &: Combes 



Maccio et al 



2005; 



20051 ) and recent theoretical work has argued that t his might even be the most 



realistic way by which galaxies get their gas (IKeres et al.l 120051 ) . Finally, in the standard 
CDM scenario, recent high-resolution cosmological simulations of galaxy formation have 
shown that angular momentum misalignment during the hierarchical struct ure formation 



can also lead to the formation of highly inclined ring/disk (IBrook et al. 



200 



Although most of the formation scenarios proposed so far seem able to account for many 
of the observed morphologies and kinematics for this class of peculiar galaxies, the relative 
weight of the two main processes, i.e. the merging of two galaxies or external gas accretion, 
remains still uncertain. The reasons for such an uncertainty mainly reside in the poor under- 
standing of many processes which are at work during gravitational interactions such as, for 
instance, gas dissipation, star formation efficiency, feedback, gas re-heating, etc., phenomena 
that all strongly affect the dynamical history of galaxies. Detailed photometric and kinemat- 
ical observations of systems in different dynamical states, both in the local universe and at 
higher redshift, can help better constraining the physics of such processes. To this aim, we 
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present a detailed study of the peculiar early-type galaxy NGC19 47, that can be considered 
as a prototype of minor-axis dust-lane galaxy (see Fig{T]). In 1978 lBertola fc Gallettal (119781 ) 
made the first study of a small sample of peculiar galaxies (including NGC1947) that all 
showed a dust-lane crossing perpedicularly the central spheroidal body: the similarities in 
observed properties let them to conclude that, by a morphological point of view, these ob- 
jects should be considered as a new class of galaxies. A larger sample of dust-lane gal axies, 



Hawarden et al 



fll98JJ), and 



including objects similar to NGC1947, have been studied by 
they found that galaxies with a minor axis dust-lane form a dominant group among all 
dust-lane galaxies. The difference of the angular momenta of the stellar, gaseous and dusty 
components is a clear proof that the three components have not all the same origin, and an 
accretion of new mat erial from outside was addressed as possible formation mechanism for 



such class of objects ( IHawarden et al 



19811 ) . 



This work is a part of an ongoing research project aimed to select a sample of peculiar 
and interacting galaxies, for which both multiwavelength photometric and kinematical data 
are available, and to compare the observed properties with the predictions derived from 
different galaxy formation scenarios. 

The main properties of NGC1947 are listed in TabO The morphological classification 
of this gala xy is somewhat controversi al. In fact, NGC1947 was initially classified as SO 



peculiar by 



de Vaucouleurs et al. 



(119761 ) and S0 3 peculiar by ISandaee fc Tammannl ( 19871 ) 



while, due to t 



r e abs ence of a photometrically detectable disk, 



Qosterloo et al. 



Bertola et al 



(jl99J) and 



(120021 ) classifie d it as an El c rossed by a series of dust lanes along the op- 



tical minor axis. As shown by 



Bertola et al. 



(119921 ) the peculiarity of NGC1947 resides in 



the presence of three parallel minor-axis dust lanes, which appear as concentric rings, in 
addition to the central one, which is typical of the dust-lane galaxies. The main dust lane 
crosses the center along the SW direction, while the others, less pronounced and parallel to 
the central one, are det ected in the NE side. The galaxy also exhibits a ring of molecular gas 
(ISage fc Gallettal Il993l ) and a disk-like distribution of ionized gas centered on the nucleus 
and at the same position angle of the dust lane. Moreover, the star s in NGC1947 rota te 
along the galaxy's major axis, perpendicular to the gas rotation axis (jBertola et al.lll992l ). 
A considerable amount of HI was detected by 



Qosterloo et al. 



)J): 



(2002): uhe HI gas is dis- 



tributed in a warped disk, which is aligned with the dust lane and it rotates in the same 
sense; the HI mass is about 3.4 x 10 8 Mq. 

This paper is structured as follows: in Section [2] we describe the observations and data 
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Table 1. General properties of NGC1947 



Parameter 



Value 



Ref. 



Morphological type 

R.A. (J2000) 

Dec. (J2000) 

Helio. Radial Velocity 

Redshift 

Distance 

Axial ratio a' 

(J ~ K) T 
Total mj (mag) 
Total rax (mag) 



SO peculiar 

05h26m47.6s 

-63d45m36s 

1100 km/s 

0.0037 

15 Mpc 

0.98 

0.931 

8.445 ± 0.02 
7.510 ± 0.03 
3.4 x 10 s 



NED^ 
NED _ 

NED 
NED 
NED 



2MASI=£ 
2MASS~ 

2MASS 




0.2 



2.65 x 10 



15.3 x 10 5 



ill 



Oosterloo et al. ( 20021 



Oosterloo et al. (2002) 



a NASA/IPAC Extragalactic Database 
b 2 Micron All Sky Survey 

reduction procedures. The morphology in the NIR and optical bands and the light and 
color distribution of the galaxy are described in Section [31 while in Section 13.11 we present 
a 2-dimensional model of the light distribution. In Section 0] we discuss the star and gas 
kinematics along the two principal photometric axes of the galaxy, and in Section [T4] is 
presented and discussed the spectral energy distribution derived for NGC1947. Finally, in 
Section [5] we draw our conclusions. 

Throughout this paper we shall adopt for NGC1947 a distance of about 15 Mpc based on 
H = 75 km s _1 Mpc -1 and an heliocentric radial velocity V = 1100 km s~ l , which implies 
that 1 arcsec = 0.07 kpc. 

2 OBSERVATIONS AND DATA REDUCTION 

Near- Infrared data - NGC1947 belongs to a selected sample of peculiar galaxies observed 
(December 2002) in the Near-Infrared (NIR) J and K bands with the SofI infrared camera 
at the ESO-NTT telescope. The field of view was 4.92 x 4.92 arcmin 2 with a pixel scale 
of 0.292 arcsec/pixel. Images were acquired in the offsetting mode: a cycle was defined by 
several images on the target, interspersed with sky frames and with an integration time of 
60 seconds; each object frame was taken with a small offset from the galaxy center and the 
sky frames were taken before and after each galaxy frame. More cycles were obtained in the 
K band than in the J band, in order to have a better estimate of the background level. A 
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total exposure time of 720 sec was obtained on the target in the J band and of 2760 sec in 
the K band. The average seeing during the observing time is about FWHM ~ 1.1 arcsec. 

The data reduction was carried out using the CCDRED package in the IRAF0 {Image 
Reduction and Analysis Facility) environment. The main strategy adopted for each data- 
set included dark subtractions, flatfielding correction, sky subtraction and rejection of bad 
pixels. Finally, all frames were registered and co-added to form the final science frames. 

(119981 ). observed at the beginning, middle and end 



Several standard stars, from 



Persson et al. 



of each observing night, were used to transform instrumental magnitudes into the standard 
J and K band systems. 

The calibrated K band image of NGC1947 is shown in Fig. [2] (right panel). 

Optical data - CCD images of NGC1947 in the Johnson V and R bands were extracted 
from the European Southern Observatory (ESO) Public Archive. They had been obtained 
with the red arm of EMMI at ESO-NTT telescope and the detector used for the observations 
consisted of two MIT/LL CCDs, with a scale of 0.332 arcsec/pixel (with a 2 x 2 binning) and 
a field of view of 22.7 x 22.7 arcmin 2 . The average seeing was FWHM ~ 0.99 ± 0.03 arcsec 
for the V band, and FWHM ~ 0.94 ± 0.03 for the R band. There were seven exposures 
available in the V band with an integration time of 100s, and five exposures in the R band 
with an integration time of 60s. The raw frames were pre-processed using the standard IRAF 
tasks for de-biasing and flat-fielding. The flat-field was performed using normalized sky flat. 
Each CCD was reduced independently: the images were dithered and shifted in order to 
obtain a complete image of both the galaxy and its surroundings. The frames of each CCDs 
were combined, after de-biasing and flat-fielding, to filter out cosmic ray hits. Finally, we 
used the IRAF task IRMOSAIC to obtain a mosaic of the images. Multiple exposures of 
identical fields were aligned using several stars present in the images and then co-added to 
produce a final median averaged image. An average value of the background emission was 
estimated in several regions of the sky far from the galaxy light. 

In order to perform the photometric calibration, several Landolt standard stars were observed 
in the same bands. 

The final, calibrated V band image of NGC1947 is shown in Fig. [2] (left panel). 
Spectroscopic data - The spectroscopic data were extracted from the European Southern 



1 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Associated Universities for 

Research in Astronomy, Inc. under cooperative agreement with the National Science Foundation. 

2 Bias frame is included in the Dark frame. 



Revisiting the formation history of the minor-axis dust lane galaxy NGC1947 7 

Observatory (ESO) public archivqfl. They have been obtained with the EMMI spectrograph 
at ESO-NTT in the RILD mode. The data, consisting of four set of spectra, were acquired 
with a slit 1.0" wide, using the GRJj5 grism with a dispersion of 5.5 nm mm" 1 , corresponding 
to 0.83 A pxl -1 , in the 380 — 700 nm wavelength range. The spectra were acquired along both 
the photometric major (P. A. = 40°) and minor (P. A. = 130°) axis of NGC1947. The total 
integration time of both spectra is 600 s; the average seeing turned out to be FWHM = 1 
arcsec. A set of spectra of standard template K0III stars were also acquired with the same 
configuration. 

Individual frames were pre-reduced using the standard IRAF image processing package; 
the wavelength calibration was made using the IRAF TWODSPECLONGSLIT package 
and a set of He-Ar-Ne lamp spectra taken for each observing night. Sky subtraction was 
performed with the IRAF BACKGROUND package, by using a narrow region at both edges 
of the slit where the galaxy contamination is minimum. Finally, all the exposures for both 
axis were co-added in the final median averaged 2D spectra, which were used to derive the 
central velocity dispersion and the rotation curve along both axes. The velo city dispersion 



and t he radial velocity w ere derived with the Fourier Quotient Technique (ISargent et al. 



1977 



Bertola et al. 



1984 ). using the standard stars as templates. One of the main source 



of statistical and systematics errors is the template mismatching; so we performed tests to 
verify if the standard star we used to derive the kinematics of the galaxy was the best one. 
These tests showed no template mismatching. In both spectra some extended emission lines 
are well visible, in particular [OIII] A5007A, H a , [Nil] AA6548,6583 and [SII] AA6717,6731. 
These lines allowed us to measure the gas rotation curves also. To measure the rotation 
curves, every record of each galaxy spectrum was fitted with a Gaussian function on the 
emission line using the task SPLOT of IRAF package. The H a line was not measured due 
to the presence of a strong absorption component, for this reason only the more extended 
[Nil] A 6583 line with a signal to noise ratio > 2 was fitted. Taking into account the average 
seeing and the EMMI pixel scale, the final stellar and gaseous rotation curves were derived by 
adopting a binning factor of 3 along the spatial direction to obtain independent data points. 
They will be presented in Fig{T2] and described in the sec]H Due to the low signal-to-noise 
ratio (which is about 20 in the center and decreases to about 10 towards the outer regions), 
the kinematic profiles are about 4 times less extended than the light profiles (Figj3j). 



A http: / / archive.eso.org/ 
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3 PHOTOMETRY: MORPHOLOGY, LIGHT AND COLOR 
DISTRIBUTION 

The most prominent feature of NGC1947 is the extended dust-lane, which crosses the galaxy 
center along the photometric minor axis. Since the dust optical depth decreases toward 
longer wavelengths, NIR photometry is necessary to reduce as much as possible the dust 
absorption, that strongly affects the starlight distribution in the galaxy. Thus, we used near- 
infrared images (J and K bands) to easily identify the inner structures of NGC1947 and 
accurate analyse the stellar light distribution. On the other hand, optical images (in the V 
and R bands) are used to derive the optical versus NIR color profiles and color maps to 
study the structure of the dust-lane. 

Isophote analysis - We used the IRAF-ELLIPSE task on the NIR images to perform 
the isophotal analysis for NGC1947 and the results are shown in FiglHJ The average surface 
brightness extends up to about 80 arcsec from the galaxy center (see Figj3j) ; in the K band, 
the half-light radius is R e = 24.4 ± 2 arcsec ~ 1.86 ± 0.15 kpc. For a semi-major axis r, in 
the range 10 < r < 40 arcsec, the ellipticity and P.A. are almost constant and equal to 0.05 
and 40°, respectively, that indicates that in this regions the isophotes are almost round and 
co-axial. For r < 20 arcsec the dust lane perturbs the ellipticity and P.A. profiles and also 
the shape parameters a^/a, which are shown in FigJHl This effect is even more evident in 
the J than in the K band. For 10 < r < 25 arcsec, the shape parameters are all consistent 
with zero, thus the isophotes do not significantly deviate from purely elliptical shape, while 
at larger distances (out to ~ 40 arcsec) they tend to be more disky (see Figj3]). In Tabj2]we 
give the integrated magnitudes and average colors within two circular apertures, derived for 
the NIR J and K bands, and optical V and R bands. The apertures were chosen in order to 
make the comparison easier with the magnitudes of NED. 

Light and color distribution - In FigJUwe show the light profiles along the major (P.A. = 
40°) and minor (P.A. = 130°) axis in the NIR and in the optical bands respectively. The 
light distributions in the optical V and R bands are very much perturbed by the presence 
of the dust-lane, in the regions between 10 and 30 arcsec on the NE side of the major axis, 
and for r < 20 arcsec on the SE side of the minor axis. Except for this, neither NIR nor 
optical profiles show any significant feature. 

We have derived the V — R (Figj5j) , J — K (FigJTl left panels) and V — K (FigJHl left panels) 
color profiles along both photometric axes of NGC1947, and the 2-dimensional V — K and 
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J — K color maps (Figj6]and FigJ7]right panels, respectively). On average, the central regions 
of the galaxy have redder colors, with a maximum value of V — R = 1.05, J — K = 1.5 and 
V — K = 5. At larger radii the colors are bluer. In all color profiles, on the NE side of the 
major axis, between 10 to 30 arcsec, double red peaks are observed, at r=18 arcsec and r=30 
arcsec where V-R=0.85 and V- K=4: such features a re due to the ring-like structure of the 



dust-lane, already suggested by 



Bertola et al. 



(119921 ) and described in SecJIJ which stands 



out very clearly in the 2-dimensional V-K color map (Figj6j). The dust-lane perturbs only 
the SE side of the minor axis color profile, within 10 arcsec from the center: this suggests 
that the dust-lane is slightly misaligned with respect to the photometric minor axis. This 
is also confirmed by the 2-dimensional J-K color map (see FigJTl right panel): with respect 
to the nucleus, which is very red, there are two "thin" and redder structures which do not 
lie on the same direction. In particular, the one on the NW side is on a parallel direction 
towards North. The full extension of the inner dust-lane is about 27 arcsec. 



Stellar population ages - We derived the integrated V-K and J-K colors in two circular 
apertures: one including the nuclear regions of the galaxy inside a radius of 1 kpc ~ 0.54i? e , 
where V — K = 3.2 and J — K = 1.02, and one out to 1.5-R e ~ 2. 85 kpc, where V 



Bruzual &: Chariot 



K = 2 .84 and J — K = 0.99. The stellar population synthesis model by 
(120031 ) were used to reproduce the integrated colors in these regions, in order to derive 
an estimate of the stellar population ages. We adopted a star formation history with an 
exponentially decreasing rate, that produces a reasonable fit of the photometric properties of 
early-type galaxies in the local Universe. It has the following analytical expression: SFR(t) = 
- exp (— t/r), where the r parameter quantifies the "time scale" when the star formation was 
most efficient. Adopting r = 1 Gyr, the correspondent evolutionary tracks were derived for 
both solar metallicity, Z = 0.02, and for a higher value, Z = 0.05; both values were assumed 
constant with age. In every model it has been assumed that stars form according to the 



Salpeterl (119551 ) IMF, in the range from 0.1 to 125M . To account for the V-K and J-K 



colors the best model is that obtained for Z = 0.05, from which we derived an age of about 
3.9 Gyrs for the inner region and of about 2.3 Gyrs for the outer one. These estimates have 
to be considered as an upper limit to the last burst of star formation. These constraints will 
turn to be fundamental in the discussion on the formation scenarios (Secj5]). 
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3.1 2-Dimensional model of the light distribution 



We performed a 2-dimensional model of the light distribution in the K band, where the 
effect of dust abso rption is weaker, and stars are accurately masked. To this aim, we used 
the GALFIT task (IPeng et al.ll2002l ). As a first attempt and according to the morphological 
classification given for NGC1947 in the previ ous work (se e Tab{T]) as an Elliptical, the galaxy 
light is modelled through a single Sersic law (jSersidll968l ): the fit is very poor, suggesting the 
presence of an additional component in the outer regions. Therefore, we adopted a double- 
compo nent model as the superposition of a Sersic law (jSersid Il968l ) plus an exponential 
profile (jFreemanlll970l ). The fit improves considerably; the fitted light profiles for the major 
and minor axis of NGC1947 are shown in FigJH]and FigfjH The structural parameters derived 
from the fit are: the bulge total magnitude m;, = 17.7 ± 0.01 mag, the bulge scale length 
r e = 74.0 ± 0.7 arcsec, the disk total magnitude rrid = 17.0 ± 0.02 mag, the disk scale length 
rh = 121.1 ± 0.3 arcsec and the Sersic exponent n = 2.83 ± 0.01, leading to a bulge total 
luminosity-to-disk total luminosity ratio, L B /L D ~ 0.43. 



Along both axes, inside a 30 arcsec radius from the galaxy center, the residuals reflect 
the presence of the dust-lane, which is more pronounced along the NE region of the major 
axis, where the galaxy appears to be less luminous than the model. Such an effect is much 
more clear in the residual image (see Fig lTOl) obtained by subtracting the 2D model to the 
galaxy. In order to analyze the structure of the dust-lane in NGC1947, we used the 2D 
model obtained in the K band, which is less perturbed by dust absorption, to match the 
light distribution in the V band. To this aim, the V band image of NGC1947 was scaled to 
the K band, by accounting for the average V-K color (see Secj3]). In the residuals (Fig JTTT) two 
main features stand out very clearly: i) the concentric ring-like structure of the dust-lane, 
superposed to a) a bright region (within 20 arcsec from the center) which is more luminous 
on the SW side with respect to the NE side. The existence of such luminous "counterpart" 
is reliable due to the fit interpolation, which is an average of the whole galaxy light (see 
Figj8]and Figj9]): this is also confirmed by the absence in the V-K color profile (see Figj6]) of 
any bluer component on the SW region, which should have appeared in the V-band residual 
image respect to the K band one. By this analysis, we estimate that the dust-lane extends 
by 49.6 arcsec from SE to NW, with a thickness of about 23 arcsec. 
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Table 2. Magnitudes and Colors for NGC1947 in circular aperture. 



aperture radius 


mj 


m K 


m v 




J-K 


V-R 


(arcsec) 


± 0.05 


± 0.04 


± 0.01 


± 0.01 


± 0.09 




15.7 


10.44 


9.35 


13.20 


11.81 


1.09 


1.39 


81.2 


8.91 


7.90 


11.26 


9.99 


1.01 


1.27 



4 STELLAR AND GASEOUS KINEMATICS ALONG THE PRINCIPAL 
PHOTOMETRIC AXES OF NGC 1947 

Major axis - The stellar kinematics is more extended on the SW side with respect to the 
NE and it is measured out to 20" ~ 1.4 kpc (Fig. [T2], bottom panel), which corresponds to a 
surface brightness of ~ 21.7 mag arcsec -2 in the V band. The rotation curve is not symmetric 
with respect to the galaxy center: in the range 2" < r < 14" NE, rotation increases, up to a 
maximum value between 55 and 80 km/s, to be compared with the corresponding distance 
range on the SW side, where almost no rotation is measured. The kinematics along NE side, 
for r > 6", is strongly perturbed by the dust-lane, as suggested by the V-K color profiles 

(FigEJ). 

The velocity dispersion remains almost constant within the uncertainties at an average value 
of ~ 133±31 km/ s, out to 10 arcsec from the center. Both the rotation velocity profile and the 



averag e velocity dispersion turns out to be consistent with those published by 



Bertola et al. 



(1992) 



The gas component too shows some rotation along this direction, and the rotation curve has 
the same extension as the stellar rotation curves (see Fig. [T2l) : it shows a very steep inner 
gradient, reaching a maximum velocity of ~ 100 km/s on the NE side and ~ 80 km/s on 
SW side, between 2" and 4" from the center. Velocities remain almost on these values up to 
the last measured p oint. These data (b oth for emission and absorption) are consistent with, 



and extend those of 



Moellenhofi ( 



1982 ). How ever, the measured rotation for the gas was not 



so evident in the data of iMoellenhofil (Il982j ). 

Minor axis - The stellar kinematics is measured out to the same distance from the 
galaxy center, i.e. 20" (Fig. [121 top panel). On average, stars do not rotate and the velocity 
dispersion has the mean weighted value of ~ 146 ± 37 km/s within 5" from the center. 
The gas rotation curve is more extended (up to about 30") and better resolved with respect 
to the stellar rotation curves. Along the major axis, the gas shows a very steep inner gradient, 
reaching a maximum velocity of ~ 200 km/s between 10" and 20" from the center. The full 
extension of the gas rotation curve is of about 60 arcsec, which is about 2 times wider than 
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the dust-lane (see Sec. [3]). 

The gas velocity curve presented in thi s work, even thoug h more extended on the SE side, 
is consistent with the one published by iMoellenhofil (119821 ) . 

The observed kinematics led us to conclude that the minor axis is the principal rotation 
plane of the gaseous disk. Moreover, for this component some rotation is observed even 
in the d irection of the major axis: t his implies that the disk is warping towards the inner 
regions (lArnaboldi h Gallettal Il993l ) and suggests that this component has not reached a 
stable configuration y et. The average velo city dispersion turns out to be consistent also 

Carollo et al.l (119931 ). according to the kinematics derived along two 



with values derived by 
intermediate directions. 

In Fig. [13] we show the stellar and gaseous rotation curves along the minor axis of 
NGC1947 superpose d on the position-velo city plot along the P. A. = 127° of the HI distri- 



bution presented by 



Oosterloo et al 



(120021 ) . The HI distribution is more extended than the 
stellar rotation curve: the HI extends out to about 90 arcsec, while the stellar rotation curve 
extends up to about 20 arcsec. The HI and gas in the dust-lane rotates in the same sense 
and with the same amplitude (see Fig flBl : this confirms th at the HI is in a warped rotating 



disk associated with the dust-lane, as already suggested by 



Oosterloo et al. 



(120021 ). 



5 DISCUSSION AND CONCLUSIONS 

In the present work we performed a detailed analysis of the structure of the minor-axis dust- 
lane galaxy NGC1947, and in what follows we will discuss the implications of our findings on 
the formation and evolution history of this peculiar object. In particular, we will discuss how 
the observed properties for NGC1947 presented in this work, compare with the predictions 
from different formation scenarios for such kind of systems. The main observed properties for 
NGC1947 that the most likely formation scenario has to account for are i) the high-inclined 
ring of gas and dust which rotates in the perpendicular direction with respect to the stars in 
the major axies of the SO-like host galaxy, and ii) the photometry and kinematics observed 
for the two decoupled components (stars and gas), in particular, the warp observed for the 
gas towards the inner regions (sec. Hj), which suggests that the inner rings of gas has not 
settled yet in the equilibrium plane; in) the estimated timescale of the formation mechanism 
has to be consistent with the age estimates of the different stellar populations. 
The presence of two components with different angular momenta, gas and dust along the 
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minor axis and stars along the major axis, suggests that NGC1947 could not be the result 
of a single protogalactic cloud collapse and that it may have experienced an interaction 
event. In this framework, the possible formation scenarios which may account for a galaxy 
with high inclined ring/disk of gas, stars and dust (discussed in the SecJT]) are i) the tidal 
accretion of gas by outside, and ii) the unequal- mass merging of two disk galaxies. In the 
following we will discriminate between these two scenarios for NGC1947. Furthermore, we 
will discuss how Smoothed Pa rticles Hydrodynamic ( SPH) simulations of isolated collapsing 
triaxial systems, presented by lMazzei fc Curirl (120031 ) . may account for peculiar galaxies like 
NGC1947. 

- Tidal accretion - 

The features observed in NGC1947 lead to the possibility of considering the galaxy as 
a tidal accretion remnant. The tidal accretion scenario accounts for a large variety of polar 
structures; many simulations show that the accretion produces warped rings and transient 



features (see 



Bournaud fc Combes 



20031 ). with processes similar to those leading to the for- 



mation of Polar Ring Galaxies (PRGs). Differently from PRGs, in the polar ring of NGC1947 
there are not detected stars associated to the gas and dust. Such Interstellar Medium (ISM) 
could be accreted by a gas-rich donor galaxy during a close passage of a pre-existing SO 
galaxy near this companion galaxy. In order to consider the accretion hypothesis, we stud- 
ied the field around NGC1947 (shown in FigJT5l) to see if there are any objects from which 
the galaxy could have accreted material: inside a radius of about 5 times the diame ter of 
the galaxy, which correspond to about 16 arcmin, (as suggested by lBrocca et al.lll997l ) there 



is the early-type spiral galaxy ESO 085-GA088, which is at the same red s hift o f NGC1947, 
with a radial velocity of 1164 km/s. As discussed by 



Qosterloo et al. 



((20021), NGC1947 

shows a warped HI disk at the same position angle as the dust lane, while the HI emis- 
sion of ESO 085-GA088 extends in a direction perpendicular to that of NGC1947: this is 
the orbital configuration ne eded to form a polar ring (of gas) through an accretion event 
( IBournaud &; Combes! 120031 ) . These observed features may suggest that the donor galaxy 
could be ESO 085-GA088: in fact, taking into account the linear relative distance of about 
16' ~ 70 kpc, and the relative velocities of the two systems, we estimate the time at which 
the interaction happened, to be about 1 Gyr ago. Such value estimated in the case of a 
parabolic orbit has to be considered a lower limit. In the case of a parabolic orbit a longer 
time is expected, which should be consistent with the upper limit to the last burst of star 
formation in NGC1947 (estimated in Sec. [3]). 
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- Merging - 

Unequal-mass merging of two disk galaxies may build a spheroidal galaxy with an high 
inclined ring/disk of gas, stars and dust. In order to test such formation mechanism, we 
will compare th e obse rved properties with those expected from simulations performed by 



Bournaud et al. 



(120051 ) . To do so, we computed the following quantities: 



the fitted light profile for this galaxy is consistent with a Sersic plus exponential lumi- 



nosity profile, wit 
• according to 



i bulge /disk=0 .43 and bulge/total=0.57 (see SecJH]); 



Bournaud et al. 



(120051 ). we evaluated the isophotal shape parameters of 
the galaxy, in the range [0.55-R 2 5th, R25th], which is the range between the bulge and the disk 
optical radius. In this range the ellipticity is almost constant, ~ 0.11, and the ~ —0.004 
suggests a bulge with a moderate boxiness; 

• we have taken into account the global kinematical properties of NGC1947 by deriving 
the v/a ratio along the major axis: in the range [0.55i?25t7i, R2Mh\ v/cr ~ 0.15, while the ratio 
Vmax/<?o is about 0.62. 



Such quantities have been compared with those p redicted by numerica 



galaxy mergers with various mass ratios (see Tab. 2 in 



Bournaud et al. 



simulations for 



20051 ): we found that 



they tu rn to be consistent w ith those expected for a merging with initial mass rati os 3:1 



or 4:1 (IBournaud et al. 



20051 ). This is a mass ratio for which 



Bournaud et al. 



2005 found 



stable and long lived polar or high-inclined rings of gas, returning to the galaxy along tidal 
tails, and equatorial rings, that appears as "dust lanes" when they are seen edge-on. Such 
unequal-mass merger has been proven to be the most likely s cenario for another peculiar 



galaxy: the double-ring ESO 474-G26 (iReshetnikov et al 



2005|) 



In the hypothesis that NGC1947 may be the result of a merging event, we compare the 
observed properties of NGC1947 with those of the galaxies that have certainly exp erienced 



Toomre 



a me rging, such as the galaxies belonging to the Toomre Sequence (hereafter TS 
19770 . 

The first step was to place NGC1947 on the TS, to determine at which merger state we 
are looking at. To this aim, we studied the position of NGC1947 on the plot of the total K 
band magnitude versus the magnitude with in a radius of 1 kpc for the galaxies belonging 



to the TS, performed by 



Rossa et al. 



(120071 ) (see Fig.25 in that paper). Thus, for NGC1947 
we estimated the total K band magnitude and the V, J and K magnitudes within a radius 
of lkpc, which are respectively Mx = —23.3 mag, my = 12.37 mag, rrij = 9.79 mag and 
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ma = 8.74 mag. The nucleus of NGC1947 lies in the region between the nuclei of the late- 
stage mergers and the normal early-type galaxies. This leads to the reasonably conclusion 
that NGC1947 could be considered a very late-stage m erger remnant. As a second step, we 



compared the V-K and J-K nuclei colors obtained by iRossa et al.l (120071 ) for the galaxies 
belonging to the TS with those derived for NGC1947 (see Secj3]), by classifying it as late- 
stage merger (FigflBland fTTI) . The J-K colors for NGC1947 are consistent with those of TS 
galaxies, while in the V-K color plot NGC1947 turns to be bluer relative to the TS galaxies: 
since no clumps of star forming galaxies are observed in NGC1947, this difference could 
reasonably due to the presence of a larger amount of dust in the TS galaxies. 
- Formation of NGC1947 through SPH simulations - 

We derive a suitabl e match of the dynami cal and photometric properties of NGC1947 



with a SPH simulation: 



Mazzei fc Curirl (120031 ) implemented SPH simulations with chemo- 



photometric evolutionary population synthesis models providing the SED from ultraviolet 
to near-IR wavelengths. The simulation includes self-gravity of gas, stars and dark matter 
(DM), radiative cooling, hydrodynamical pressure, shock heating, artifici al viscosity, star 



Mazzei fc Curir 



forma tion and feedback from evolved stars and type II SNe, as described in 
(120031 ). The starting point is a triaxial collapsing system initially composed of 30000 parti 
cles, 15000 of gas and 15000 of DM with a relative mass ratio of 0.01. Th e total mass, M, 



Mazzei fc Curir 



of the collapsing system is 1O 13 M0. The system is built up as described by 
((20031), i.e., with a spin parameter, A, given by | J| \E\°' 5 /(GM ' 5 ), where E is the total en- 
ergy, J the total angular momentum and G the gravitational constant, equal to 0.06, and 
it is aligned with the shorter principal axis of the DM halo. It is also assumed a triaxiality 



ratio of the DM halo, r 



b 2 )/(a 2 — c 2 ) = 0.84, where a > b > c, whic h corresponds to 



Mazzei fc Curir 



an ave rage radius of about 1 Mpc. Under these conditions, as found out by 
((20031), the morphology of the resulting galaxy is highly related to the initial properties 
of the halo, which drive the galaxy formation and evolution. In particular, the dynamical 
timescales of stars and gas are quite different and, as a consequence, the two components 
become decoupled during the evolution of the system. 

Here, we present the synthetic SED which accounts for chemical evolution, stellar emis- 
sion, internal extinction and re-emission by dust in a self-consistent way, as described in 



((Mazzei et al. 



1992 



1995( 1. This allows us to extend the SED over four orders of magnitude 



in wavelength, i.e., from 0.1 to 1000 /im. So the model self-consistently provides morpholog- 
ical, dynamic and chemo-photometric evolution. Fig. [T3] compares the predicted SED with 
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the available data for NGC1947. Such a model suggests that star formation started about 
14.5 Gyrs ago, and that the maximum was reached about 9 Gyrs ago. The average age for 
the inner stellar generations is dated at about 3 Gyrs: t his could be due to the still efficient 



inflow of the gas and is consistent with the findings by ISerra et al.l (120071 ) , which included 
such a galaxy in the sample of centrally-rejuvenated objects. 

Dynamical predictions of such SPH simulations agree with kinematical features observed 
for NGC1947 (see Sec. H]): in particular, the model predicts a maximum rotation velocity of 
about 20 km/s for the stars and of 120 km/s for the gas, and a velocity dispersion in the 
central regions of about 130 km/s. 

- Conclusions - 

Putting all the above evidences together, we now try to address which are the observa- 
tional aspects that can help to disentangle in a non ambiguous way the possible formation 
scenarios for NGC1947 discussed above. 

We have found that an unequal-mass merging of two disk galaxies, with a mass ratio 
in the range 3:1 - 4:1, could account for the average structure and kinematics observed for 
NGC1947. But this scenario fails to reconcile i) the epoch of the merging event with the 
estimates of stel lar population ages in NGC1947; and ii) the gas kinematics. According to 



the simulations (BournaudetaL 



20051 ). after 3 Gyrs from the merging event, the outer re- 
gions of the remnant are characterized by typical signs of the past interaction, like tidal tails 
and/or shells. Such features are not observed in NGC1947. To account for this, the merger 
would have occurred very long ago (~ 10 Gyrs) and the ring should be very long-lived: 
these predictions are not consistent with both i) the upper limit to the epoch at which the 
interaction happened for NGC1947, i.e. ~ 3 Gyrs (see Secj3]) and ii) the warped structure 
of the ring, which suggests that the polar component has not reached a stable configuration 
yet (see secH}. Therefore, we can reasonably rule out the merging scenario for the formation 
of NGC1947 and we can regard the recent polar accretion from a gas rich donor galaxy as 
one of the most probable interaction event in the formation history of this galaxy. 
We have also discussed how such a peculiar object may be the result of processes of galaxy 
formation and evolution under particular initi al conditions, linked to the ratio between bary- 



onic and dark matter: the SPH simulations by lMazzei fc Curirl (120031 ) can reasonably account 



for the observed photometry and kinematics of the galaxy. We aim to further investigate 
in detail such aspect by comparing the predictions of such a model and observations, but 
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Figure 1. Color composite image of NGC1947 assembled from EMMI-NTT images in the V (blu channel), J (green channel) 
and K (red channel) bands. The North is up, while the East is on the left of the image. 

they are beyond the aims of the present work and they will be the subject of a forthcoming 
paper. 
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Figure 10. 2D fit of NGC1947. Left panel - K band image of NGC1947. Middle panel - Model of the galaxy. Right panel - 
Residual image. 




Figure 11. Residuals of the subtraction of the K band model to the V band image. 
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Figure 12. Stellar (filled points) and gaseous (triangles) rotation curves along the major (bottom panel) and minor axis (top 
panel) of NGC1947. The dashed line is the velocity of the galaxy center. In the bottom-right corner is plotted the error bar 
that represents the mean error on the velocity of the gas. 
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Figure 13. Stellar (filled points) and gaseous (triangles) rotation c urves along the minor axis of NGC1947 superposed on the 
position-velocity plot of the HI along the P. A. = 127o presented bv lOosterloo et al.l (2002). 
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Figure 14. Continuous line shows the prediction of our model (see text); filled squares are V, J, and K integrated data given 
in sec. [3] the othe r data are from NED, in particular IRAS data at 12, 25, 60 100 fim and ISOPHOT data at 170 ^m of 
IStickel et all (|2004h . 
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Figure 15. Image of a field of 25 arcmin around NGC19 47. The two boxes on the left show the HI distribution for NGC1947 
and ESO 085-GA088 published bv lOosterloo et alj <2002fi . 
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Figure 16. V-K color measured within lkpc apertures as a function of the Toomre sequence merger stage. The dotted horizontal 
line indicates the V-K color for a SSP W 10 yr population of solar metallicity, which provides a reasonable fit to the measured 
colors for early-type galaxies. 
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Figure 17. J-K color measured within lkpc apertures as a function of the Toomre sequence merger stage. The dotted horizontal 
line indicates the J-K color for a SSP 10 10 j/r population of solar metallicity. 
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